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Abstract 

The deconfinement transition region between hadronic matter and quark- 
gluon plasma is studied for finite volumes. Assuming simple model equations 
of state and a first order phase transition, we find that fluctuations in finite 
volumes hinder a sharp separation between the two phases around the critical 
temperature, leading to a rounding of the phase transition. For reaction 
volumes expected in heavy ion experiments, the softening of the equation of 
state is reduced considerably. This is especially true when the requirement of 
exact color-singletness is included in the QGP equation of state. 
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I. MOTIVATION 



A primary goal of relativistic nuclear collisions is the observation of a phase transition 
of confined, hadronic matter to a deconfined quark-gluon plasma. One of the proposed 
signatures, namely hydrodynamic flow, is based on the presumed softening of the eqution 
of state due to the rapid increase of the entropy density. It has been investigated in the 
framework of relativistic hydrodynamic models The expansion of once compressed 

matter is predicted to be delayed in the case of a QGP, which in turn leads to a reduction of 
the transverse (directed) flow [|. This is mainly due to the fact that the sound velocity 
vanishes for energy densities in the mixed phase. A smooth crossover transition within an 
assumed interval of AT = O.lTc, on the other hand, results in drastically reduced time 
delays as compared to a sharp transition [^. 

As is well known, a rounding of sharp first order phase transitions is expected due to 
explicit finite size effects 0. The importance of fluctuations of the coexisting phases in 
small volumes of strongly interacting matter has already been pointed out in [§] for the 
case of a liquid-gas phase transition, and in for the deconflnement phase transition. As 
a consequence, it was claimed that the observation of two separate phases in heavy ion 
collisions might be hindered. In this work we explore this behaviour in more detail, starting 
from rather simple model equations of state. We put special emphasis on the question, how 
the requirement of color singletness of the QGP phase affects the phase transition on top of 
the fluctuation effect. 

The limited reaction size of a heavy ion collision is a generally ignored problem of the 
experimental search for a quark-gluon plasma. According to one-fluid dynamical model 

^However, in a three-fluid hydro dynamical model |^ the directed nucleon flow is already lowered as 
compared to the usual one-fluid models (which assume instantaneous local thermalization between 
projectile and target). It will be exciting to learn whether the softening of the EOS is also signaled 
in this model. 
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calculations, one expects to produce deconfined matter in volumes of about 10 — 50 fm^ 
at fixed target energies p[. Similar values of typical QGP cluster sizes are found within 



a microscopic approach in combination with a percolation model |[T^. At collider energies 
the relevant reaction volumes for heavy systems are estimated to be much larger: assuming 
an initial interaction volume of Vq = 7r(1.15A^/^)^ro with a formation time Tq ^ Ifm/c and 
subsequent expansion one may expect (at the onset of confinement) V ~ 750 fm'^ for Au+Au 
collisions at RHIC energies {^/s = 200 AGeV) and V ^ 1350 fm^ for Pb+Pb at LHC energies 
(v/i^5.5ATeV) 111. 

However, the longitudinal flow velocities exceed the thermal motion by far. Thus, only 
subsystems of smaller volume can be regarded as being in approximate "global" thermal 
equilibrium. Only the latter are suitable for the study of finite size effects. Therefore, if we 
require the local thermal motion to be of the same order as the relative flow velocity, one 
must restrict the study of finite size effects to regions of about one unit of rapidity (cf. [0)- 
The volume of such (possibly multiple) subsystems would thus only be < 100 fm^ for both 
RHIC and LHC. 

Another motivation for the study of finite size effects are lattice QCD calculations, which 
have to be extrapolated from a finite lattice (although with periodic boundary conditions) 
to infinite matter |1T3| , |T^ . Depending on the number of grid points in the spatial direction, 
the presently accesssible volumes are in the range of ~ 50 — 100 fm'^. The extrapolation to 



infinity is not clear Our investigation can yield some insight what can be phenomeno- 
logically expected when going from finite volumes to the infinite volume limit. 



II. PHASE COEXISTENCE AND FLUCTUATIONS 

Let us consider a finite volume of strongly interacting matter in the grand canonical 
ensemble with vanishing net-baryon density. The model equation of state is constructed by 
matching a low energy density phase (a hadron gas) and a high energy density phase (the 
quark-gluon plasma) : For infinite matter the system undergoes a first order phase transition 
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at a critical temperature Tq- However, when we assume a finite volume V of the system, 
statistical (thermodynamical) fluctuations are not negligible; they do not allow for a sharp 
transition between the two phases in this case. In general, the probability p of finding the 
system in a state x is given by p(x) ~ exp (— /5F(x)), where F{x) is the free energy of 
the system [^,^ (see also the appendix). Let us introduce an "order" parameter ^ for the 
quantitative characterization of the macroscopic state of the total system, such that ^ = 1 
corresponds to the pure hadron phase and ^ = corresponds to the pure quark-gluon phase. 
For the mixed phase the fractional volumes are defined as Vh = ^V^ and Vq = (1 — ^)V. 

It is now assumed that the partition function of the total system factorizes into the 
partition functions of the two individual phases for fixed ^ (see appendix). Then there are 
no correlations between the possible microscopic energy eigenstates of the two subsystems 
(the two phases). The free energy of the total system can thus be expressed in terms of 
the free energy densities of the individual phases /q and /h and the volume fraction of the 
hadron phase ^: 

F^{T, V) = [A(T, ^V)^ + /q(T, (1 - 0V){1 - m (1) 
The normalized probability density is a function of the order parameter: 

^^^^ /oMeexp(-/3F^(T,\/)) 
The expectation value of any intensive thermodynamic quantity, A{T, V), is then calculated 
as: 

A{T,V)= f'd^pi^;T,V)[Aj,iT,^V)^ + A,{T,il-OV)il-0] . (3) 

A reminder of a thermodynamic derivation of (^ and (§) can be found in the appendix. 

As long as explicit finite size effects in the equations of state of the individual phases are 
ignored (i.e. the free energies Fj are linear functions of Vi), the free energy densities equal 
the negative pressure within the two phases: 
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When inserting this into (|l]) and (Q), respectively, one sees immediately that for large vol- 
umes, V, of the total system and a given temperature, T, the occurence of the phase with 
the lower pressure is strongly suppressed. For V oo the coexistence of both phases re- 
quires Ph = Pq, which is the Gibbs condition for mechanical two phase equilibrium. This 
holds also for the case of explicit finite size modifications of the equation of state, because 
the correction terms vanish by definition in the limit V ^ oo. 



For illustration let us consider first the case of infinite matter equations of state, i. e. 
without explicit finite size corrections. The constituents of the low energy density phase 
are well-established non-strange hadrons up to masses of 2 GeV. The system is treated 
as a mixture of relativistic, non-interacting Bose-Einstein and Fermi-Dirac gases. As the 
net-baryon number is zero, the free energy reads 



where the + stands for fermions, the — for bosons and gi denotes the spin and isospin degen- 
eracy of particle species i. To take into account repulsive interactions, all thermodynamic 
quantities are corrected by the Hagedorn factor 1/(1 -|- e/AB) [|TB|, where e is the energy 
density of the point particles and B is the bag constant. The bag constant is chosen as 
B = 200 MeV/fm^ for all calculations presented here. 

For the equation of state of the high energy density phase (the quark-gluon plasma) we 
take a simple bag model EOS for massless quarks, their antiquarks, and gluons in an MIT 
bag of infinite volume, where the number of quark flavors is A^^q = 2: 



III. THE EQUATIONS OF STATE FOR INFINITE MATTER 






(6) 




{A,{T) - MT)) + MT) . (7) 
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Note that fh{Tc) = fq{Tc). The expectation value of an intensive thermodynamic quantity 
hke pressure, entropy density or energy density at a temperature T ^ Tq can be calculated 
with this formula. 

Fig. |I] (top) shows the expectation value of the order parameter < ^ >, the average 
hadronic volume fraction, as a function of the temperature for different system sizes. This 
quantity is calculated as 

<e>(T,r)= f'd^pi^;T,V)^ . (8) 

For volumes V = 100 fm'^ or less, < ^ > (T) deviates clearly from a simple step function. 
Even below Tq, the system is not purely hadronic. The quark phase has a finite probability 



(see also [[T^])- On the other hand, the hadronic phase contributes, with similar probability, 
even above Tq. Note that the order parameter in our model is simply the volume fraction 
of the hadronic phase. Thus, for a volume of 25 fm^ and a temperature of 10 MeV above 
Tq one can read off Fig. ^ that the system is composed of about 10 % hadron gas and 90 % 
QGP (this admixture for finite sizes is different from the predicted persistence of clusters in 



an infinite plasma [jT^)- The unfavorable phase is suppressed, according to (0), due to the 
differences of the free energy. However, any composition of the total system (any value of 
^) must correspond to a finite total free energy content, because the volume is finite. Thus, 
the suppression relative to the state with minimum free energy content is finite, i. e. the 
probability is non-zero. 

Fig. |l| also shows the energy density e/T^ and the entropy density s/T^ as a function 
of temperature for different volumes of the system. In the case of an infinite volume, the 
first order phase transition is reflected by the sharp discontinuity of both quantities at 
Tq = 143 MeV. Due to the finite probability of fluctuations, the average values of the 
energy density and the entropy density at a temperature T are different for finite systems. 
Below Tc, the presence of the quark phase, although a small fraction of the total volume, 
increases the energy density and the entropy density. Above Tc, the contribution of the 
hadronic phase still lowers the average values of e and s leading to the expected rounding of 



the phase transition. A functional form of thermodynamic quantities has been parametrized 
PJT9[] ; this was used to model the assumed smooth crossover transition and for studying the 
resulting physics. In any case, for finite size systems such a crossover has to be expected. 
Here we calculate such a behaviour for strongly interacting, finite systems on the basis of a 
very simple model without free parameters. 

Fig. ^ now shows the (hydrodynamically relevant) ratio P(e)/e vs. e for various system 
sizes. Here the 'softest point' of the equation of state vs. V are the respective minima of 
the curves. The thermodynamics of the system exhibits a very distinct volume dependence 
with respect to this quantity. Observe the strong influence of V over a wide range of energy 
densities. This could not be clearly seen in Fig. ||, where the transition appears within a 
rather narrow region of temperature. 

It has been suggested that a clear peak of the lifetime of the mixed phase as a 
function of the collision energy could signal the QCD phase transition in heavy ion reactions. 
This peak should be observable in a particular window of the initial energy density, around 
the softest point of the equation of state. The comparatively small pressure prevents a fast 
expansion of the system around that well located point. 

However, as can be seen in Fig. 0, the minimum of P(e)/e is much less pronounced for 
finite volumes than for an infinite volume. At low temperatures — which correspond to 
low energy densities — the influence of a small admixture of quark matter on P and e is 
particularly strong, because the absolute values of P and e are very large in the quark phase 
in the simple MIT bag description {—P = e = -B for T = 0) as compared to the values in the 
hadron phase (P = e = for T = 0). In Fig. |]full circles on the curves indicate the values 
of energy density, where the pressure of the pure quark phase would become zero according 
to the simple bag model equation of state, eq. (^. Whether the concept of thermodynamic 
fluctuation theory is adequate for lower energy densities and how the equation of state for 
the supercooled quark matter phase might look like at temperatures considerably lower than 
Tq remains to be investigated. 
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IV. QUARK MATTER EQUATION OF STATE WITH COLOR SINGLET 

CONSTRAINT 



A grand canonical partition function for a quark-gluon plasma droplet with the necessary 



requirement of color-singletness was proposed in |20|. The internal SU(3)c-symmetry is 



accounted for by applying a group-theoretical projection technique on color singlet 
states to the grand canonical partition function of a noninteracting quantum gas of massless 
quarks and gluons. The authors find that, due to the color confinement, the internal degrees 
of freedom in a finite plasma droplet are effectively reduced as compared to the infinite 
matter equation of state. Also, the finite level density (in a finite system) for the single- 
particle eigenstates lead to an additional effective reduction. This is explicitely accounted 
for. In the following, we use the corresponding free energy as the equation of state for the 
high density phase. It reads ||20| 



F(T, V) = BV -T {X-Y) + '^-T {\og{D) - log(7r)) + T (4 log(C) + log(2 v^vr)) (9) 



\60 45; 



Y = ^ [\n^+^^] RT 



9 ^ 27/ V30 ^45/ 

where R= (^)^''^ is the radius of the spherical plasma drop and Nq = 2 denotes the number 

of (massless) quarks. This equation of state (for the quark-gluon plasma) incorporates an 

dF F 

explicit volume dependence of the pressure at fixed T. In this case, -Pq(T) = —^77- 7^ — tj^- 
Therefore, eqs. @) and (^) have to be evaluated numerically. The partition function is 
derived in saddle-point approximation which breaks down for RT 0. It agrees within 



30% with the exact value for RT ^ 1 |2^. Therefore, we approximate the free energy 



density by f{R) = /(1/T) for plasma volumes (1 — C,)Vtot, which correspond to spherical 
droplets of sizes less than R < 1/T. The hadronic equation of state of Sec. |T| is used, 
without any explicit finite size modifications. 

Fig. ^ (top) shows the free energy density as a function of the order parameter ^ (the 
hadron volume fraction) using the color singlet equation of state for the quark phase. The 
temperature is chosen to be T = Tc, where we define the critical temperature at the point 
where e(Tc) = (eg(Tc) + e/i(Tc))/2. Thus, Tc denotes the temperature, for which both 



phases contribute with equal probability to the state of the total system. In Sec. |T| the 
critical temperature was found to be independent of the volume — only the width of the 
transition was affected by the flucutations. 

For the finite size corrected equation of state, however, the critical temperature shifts 
to higher values at finite volumes (see below). As one expects, the finite size modifications 
of the equation of state affect the free energy density most strongly for small systems. The 
infinite volume limit converges to a constant free energy density. Fig. |^ (bottom) shows 
the resulting probability densities as functions of the order parameter ^. A 'two-hump' 
structure is observed, thus favoring the dominance of one of the phases against a two phase 
composition with equal volume fractions for both phases. Such a two-hump structure is 
generally expected to occur for first order phase transitions 0. 

One should emphasize that the peaks at ^ = and ^ = 1 are most pronounced for 
the largest volume. This is because the free energy, and not the free energy density, enters 
eq. (H) for the calculation of the probability density. The physics of large systems is therefore 
extremely sensitive to small variations of the free energy density. Because of this very distinct 
hump structure, either one phase or the other is present at Tc (which results from the explicit 
finite size corrections). Then, a first order phase transition occurs dynamically in a way that 
the high temperature phase (or the low temperature phase, respectively) is supercooled (or 
superheated, respectively), up to a point where bubble nucleation starts to convert the now 
unstable phase to the more stable phase . 



Fig. ^ shows the temperature dependence of the order parameter, the energy density and 
the entropy density, as in Fig. |1], but now for the quark matter equation of state with the 
color singlet constraint. The main difference is the shift of the critical temperature to higher 
values for smaller volumes. For systems of finite size the modified equation of state yields 
a lower pressure at a given temperature, because the internal degrees of freedom are grad- 
ually 'frozen' with decreasing volume ||20|. Thus, even when flucuations are neglected, the 



mechanical Gibbs equilibrium between the two pressures would be reached for temperatures 
Tc > Tq'. Here Tq' stands for the critical temperature of the infinite system. 

Important consequences follow for the bulk quantities e/T^ and s/T^: the latent heat 
and the jump in the entropy density are considerably reduced for small systems. First, the 
effective number of degrees of freedom in the hadron phase increases with temperature: the 
restmasses become more and more negligible as compared to the kinetic energies. Secondly, 
as mentioned before, the effective number of degrees of freedom in the quark phase is reduced 
for finite volumes due to the requirement of color singletness. As stated in |T^, the active 



degrees of freedom, which are "quantified" as s(T), determine completely the gross behaviour 
of the thermodynamics near the phase transition. Fig. ^ also shows that the smearing 
due to the fluctuations is less pronounced than for the infinite matter equation of state in 



Sec. |T|. The rounding effect of a fluctuating phase composition of the system appears to 
be counteracted by the effect of the volume dependent reduction of quark-gluon degrees of 
freedom. 

Fig. 1^ depicts the volume dependence of the shift of the critical temperature, ATq = 
Tq — T^. The temperature shift exhibits an approximate power law, ATc ~ l/V^ with 
a ~ 0.7. The analysis of Binder and Landau 0, who examined the first order phase 
transition of a finite system in the Ising ferromagnet model, advocates a coefficient a = 1. 

As in Fig. ^ Fig. |^ shows the ratio P(e)/e vs. e within the present scenario. Again, 
those values of the energy density are marked, which correspond to pressure zero of the 
pure quark phase. In clear contrast to Fig. 0, the 'softest point' is now characterized by less 
pronounced minima of P(e)/e for reasonably small system sizes. For volumes of < 25 fm^, 
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the minimum of the curve vanishes completely. Hence, the lifetime signal |l]J^, which is based 
on hydrodynamic considerations and infinite matter equations of state, will be extremely 

damped in a more realistic scenario of heavy ion collisions. 

dp 

Fig. 1^ shows the speed of sound (squared) function of the energy density 

for three different cases. For infinite volume, the phase transition is truely first order. This 
can be seen from the vanishing speed of sound in the mixed phase. This is the cause of the 
pronounced time delay in hydrodynamical simulations: even if there are strong gradients in 
the energy density, the mixed phase cannot perform mechanical work. Therefore, it does 
not expand on its own account. Deflagration fronts with small velocities convert the mixed 
phase into hadrons, leading to slow cooling and expansion p. 

A finite system of volume V = 100 fm'^ is also depicted in Fig. |^ for the two different 
equations of state of the QGP. Neglecting the color singlet constraint (but including the 
effect of thermodynamic fiucutations in finite volumes) leads to a strong reduction of the 
speed of sound in the transition region. Thus, one still might expect a moderatively long 
lived fireball of the mixed phase. However, the speed of sound does not vanish in a sharp 
region of energy density. Rather, it is smoothly damped in the large e domain between 
20MeV/fm^ and 200MeV/fm^. This effect would vastly smear out the lifetime signal if one 
scans the flow excitation function over a wide range of collision energies. 

After the requirement of color singletness is included, cj. is much less reduced. The dip is 
shifted to higher (factor ~ 2) energy densities. Although the rather sharp breakdown of the 
speed of sound in a well located energy density domain is recovered, the hydrodynamical 
expansion solutions will now look much different from the infinite matter scenario: The 
rather high values of must lead to a much more rapid expansion than for infinite matter 
equations of state. The implementation of the present results, based on finite volumes 
and global equilibrium and also supercooling, into hydrodynamical calculations, remains a 
difficult enterprise for future investigations. 
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V. CONCLUSION 



We have discussed the thermodynamic bulk properties for finite systems at the phase 
transition and the detailed behaviour of the free energy density as a function of the tempera- 
ture. For finite volumes, V < 100 fm'^, corresponding to the expected plasma volumes, there 
is a considerable rounding in the variables e/T^ and s/T^ around Tc- As a consequence, 
the speed of sound does not vanish in the mixed phase. This is infered, under simple as- 
sumptions, from basic thermodynamic considerations. Fluctuations of the two phases in a 
finite system lead to a smooth transition between the low temperature regime — where the 
hadronic phase dominates the system — and the high temperature regime — where the pure 
quark phase is most probable. 

The requirement of exact color-singletness within the quark phase leads to a shift of the 
critical temperature to higher temperatures for finite volumes. We observe a double hump 
structure in the probability distribution for the actual phase composition during the phase 
transition. In the model scenario the speed of sound is considerably increased in the mixed 
phase, if the requirement of color-singletness is taken into account. The significance of the 
time delay signal |jl],0 for the experimental detection of a QGP phase in heavy ion becomes 
questionable. 

In addition, our investigations also show intuitively that the extrapolation from the be- 
haviour of a finite system to its infinite volume limit may, in fact, be rather model dependent. 
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APPENDIX A: THERMODYNAMIC DERIVATION OF THE MODEL 



For a canonical ensemble all thermodynamic quantities like the free energy F, the entropy 
S, the pressure P and the energy E are determined by the partition function Z: 

Z = Tr{e-^^} = ^e~^^" , F = -rinZ . (Al) 

n 

Once F{T, V) is fixed, S, P, E follow from standard relations. 

The main assumption for the case of two coexisting phases is separability of the energy 
spectra E"^*^ of the single phases i. It is assumed that the energy eigenvalues of the total 
system E^ factorize: 

m) = ES{0 + E^nHO ■ (A2) 

Here ^ is a not yet specified order paramter which characterizes the configuration of the 
total system, i. e. the relative importance of the individual phases. The partition function 
for a given ^ of the two phase system then reads: 

ni n2 

F(0 = -T In Z(0=F« (0 + ^(2^0 • (A3) 



(^3) guarantees that the simple factorization in ([A3[) holds. Since ^ is a (continuous) order 



parameter, any value of ^ corresponds to a different state of the two phase system. As in our 
intuitive choice for in the main text, we now assume that, without loss of generality, ^ can 
take any number in the range between and 1. The total partition function of the system 
for a discretized spectrum of the order parameter and a given total volume V^°^ reads: 

z^z(ei) + z(6) + ...z(e„) ,with 6 = ^^^ • (A4) 

n 

The probability for the system being in the state is then given by 

Ap(e.) = ^ , EAp(6) = l . (A5) 
^ i=i 
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Prom this we infer the probabihty density of the continuous case: 

p(^) = ^ii^_ . (A6) 

We now evaluate the free energy of the total system as well as the other thermodynamic 
quantities 5"*°*, P*°* and E^^. In the discretized case we have 

n n 
i=l i=l 

Differentiating with regard to 13 (as an independent paramter) yields: 

F''" = T.^^^r—mi^m) LpiOFmc . (A8) 

The entropy is readily obtained as 







^tot 



V 



\v 



rp n 



In z + i V ^— ^Frn - V I 



^ r^^^^ ("^) '^^^ ^ I'piOSmC , (A9) 
similar for the pressure, 



ptot 



T — 77 TT \T 



V J ' z V y ^ 



/^(o (-^) iTde = J\{Opm^ , (Aio) 



and the energy 
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E 



tot 



jptot _j_ rptot ^tot 

n 



i=l 

r-1 







p(0 



F(0-T 









JO 



(All) 



)- (|All| ) do correspond to the prescription eq. stated within the main text. 
We now choose the order parameter ^ as the volume fraction of one of the two phases. 
The free energy of the system as a function of the order paramter thus becomes (confer to 
eq. im 



F{o = \f^'\^r"')^ + /(2)((i - ovm^ - 



V 



tot 



(A12) 



where is the free energy density of phase i. Note that the individual free energy densities 
can depend on the volume explicitly. For infinite matter equations of state, however, this is 
not the case. We now get for the other quantities: 



V 



sW(^V'°')^ + s(2)((l - - 

where s*^*^ denotes the entropy densities, 

p(i)^^^tot^^ + P^^\ii - - 



tot 



(A13) 



\ dV J 



(A14) 



and 



EiO = F{0 + TiOSiO = [e^'\^V'°')^ + e(2)((l - OV^i^ - 
where e*^*-* are the energy densities. 



V 



tot 



(A15) 
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FIG. 1. The order parameter (expectation value of the volume fraction of the hadron gas 

< ^ >) vs. the temperature for systems of different sizes (top). The energy density e, divided 

by the temperature T to the power of four (middle), and the entropy density s, divided by the 

temperature to the power of three (bottom), are also shown. The bag constant is B^^^ = 200 MeV. 

Equations of state for infinite matter (hadronic and QGP) are used. 
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FIG. 2. Pressure P divided by the energy density e vs. the energy density for systems of 
different sizes. The bag constant is 5^/^ = 200 MeV. Equations of state for infinite matter 
(hadronic and QGP) are used. 
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FIG. 3. Free energy density / as a function of the hadron fraction for systems of different 
sizes at Tq (top) and the resulting probabihty densities p(^) (bottom). The bag constant is 
SV4 = 200 MeV. The color singlet constraint is taken into account for the QGP equation of 
state. 
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FIG. 4. The order parameter (expectation value of the volume fraction of the hadron gas 
< ^ >) vs. the temperature for systems of different sizes (top). The energy density e, divided 
by the temperature T to the power of four (middle), and the entropy density s, divided by the 
temperature to the power of three (bottom), are also shown. The bag constant is B^l^ = 200 MeV. 
The color singlet constraint is taken into account for the QGP equation of state. 
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FIG. 5. Shift of the critical temperature ATq vs. the systems size V. The bag constant is 
= 200 MeV. The color singlet constraint is taken into account for the QGP equation of state. 
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FIG. 6. Pressure P divided by the energy density e vs. the the energy density for systems of 
different sizes. The bag constant is B^l'^ = 200 MeV. The color singlet constraint is taken into 
account for the QGP equation of state. 
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FIG. 7. Speed of sound (squared) as a function of energy density e for three different cases: 

1) infinite volume of the system (full line). 

2) V = 100 fm^ using the infinite matter EOS (dotted). 

3) V = 100 fm^ using the EOS with color singlet constraint (dashed). 
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